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We investigate the nonhnear dynamics of inertia-gravity (IG) wave modes in three- 
dimensional (3D) rotating stratified fiuids. Starting from the rotating Boussinesq equations, 
we derive a reduced partial differential equation, the GGG model, consisting of only wave- 
mode interactions. We note that this subsystem conserves energy and is not restricted to 
resonant wave-mode interactions. In principle, comparing this model to the full rotating 
Boussinesq system allows us to gauge the importance of wave- vortical- wave vs. wave- wave- 
wave interactions in determining the transfer and distribution of wave-mode energy. As in 
many atmosphere-ocean phenomena we work in a skewed aspect ratio domain H/ L {H and 
L are the vertical and horizontal lengths) with Fr = Ro < 1 such that Bu = 1, where 
Fr, Ro and Bu are the Froude, Rossby and Burger numbers, respectively. Our focus is on 
the equilibration of wave-mode energy and its spectral scaling under the infiuence of random 
large-scale (fc/) forcing. We present results from two sets of parameters: (i) Fr = Ro w 0.05, 
H/L = 1/5, and (ii) Fr = Ro « 0.1, H/L — 1/3. As anticipated from prior work, when 
forcing is applied to all modes with equal weight, with Fr = Ro ~ 0.05 and H/L = 1/5, the 
wave-mode energy of the full system equilibrates and its spectrum scales as a power-law that 
lies between and k~^^^ for kf < k < k^, where fc^ is the dissipation scale. For the same 
parameters, when forcing is restricted to only wave modes, the wave-mode energy fails to 
equilibrate in both the full system as well as the GGG subsystem at the resolutions we can 
achieve. This cleary demonstrates the importance of the vortical mode (by facilitating wave- 
vortical-wave interactions) in determining the wave-mode energy in the rotating Boussinesq 
system. Proceeding to the second set of simulations, i.e. for the larger Fr = Ro ^ 0.1 in a 
less skewed aspect ratio domain with H/L = 1/3, we observe that the energy of the GGG 
subsystem equilibrates and is resolution independent. Further, the full system with forcing 
restricted to wave modes also equilbrates and both yield identical energy spectra. Thus it is 
clear that the wave-wave-wave interactions play a role in the overall dynamics at moderate 
Ro, Fr and aspect ratios. Apart from theoretical concerns addressed in the Conclusion, from 
a practical standpoint these results highlight the difficulty in properly resolving wave-mode 
interactions when simulating realistic geophysical phenomena. 
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I. INTRODUCTION 

Inertia gravity (IG) waves, resulting from the rotating and stratified nature of geophysical fluids, 
play an important role in atmosphere-ocean dynamics [Ij. A broad overview of their properties 
with relevance to the ocean can be found in Garrett & Munk , and a recent review focussing on 
the observational characterization of the oceanic wave- field can be found in Polzin and Lvov [3]. 
The general dynamics of these waves are reviewed by Sommeria & Staquet [3]. Further, Wunsch 
& Ferrari [5j place these waves in context when considering the general circulation of the ocean, 
while Fritts & Alexander [6] review their influence on diverse phenomena in the middle atmosphere. 

In their influential papers, Garrett & Munk [2], [7] pointed out the importance of a deeper 
understanding of the dynamics of the IG wave modes. In this regard, one avenue of progress 
has been to focus upon the subset of resonant interactions among these waves [8j. In general, 
the dispersive nature of the linear waves results in reduced nonlinear transfer between modes by 
dispersive phase scrambling. However, resonant interactions are special nonlinear interactions 
for which phase scrambling is absent, and nonlinear transfer remains strong. Thus resonant 
interactions are critical for the dynamics when they are present. See for example |9j for a broad 
introduction to nonlinear dispersive equations, [TTT for statistical theories of the so-called weak 
turbulence, [H], [l2] for geophysical perspectives and [I3],[I3] for the roles of resonances and 
near-resonances in the specific case of the 3D rotating Boussinesq system. 

Investigating resonant IG waves, McComas & Bretherton [8] distinguished between three 
classes of resonant interactions that lead to induced diffusion, elastic scattering and parametric 
sub-harmonic instability. The picture put forth was, for a large-scale initial data: (i) elastic 
scattering leads to 3D isotropy of an initial condition that is only horizontally isotropic, (ii) 
sub-harmonic instability moves energy downscale, (iii) induced diffusion tends to force the high 
wavenumber spectrum towards equilibrium. A detailed account of these resonances can be 
found in the extensive review by MuUer et al. [15]; their review also describes approaches that 
account for more than resonant interactions, but involve other approximations (such as the direct 
interaction approximation). Further, the use of resonances to derive kinetic equations under the 
weak turbulence paradigm has been an active area of work — for an overview see Lvov et al. |16] . 
In particular recent work has shown a correspondence between observed IG energy spectra and 
particular steady state solutions of the relevant kinetic equations |17j.|18]. 
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Unfortunately, as was pointed out by McComas & Bretherton [8], resonances do not account 
for the bulk of interactions among IG waves. A formal statement regarding the sparsity of these 
resonances can be found in Babin et al. [T9] . This is all the more an important issue when i?o, Fr 
are small (implying strong rotation and stratification) but far from zero, i.e. ^ i?o, Fr < 1 
(where Ro,Fr are the Rossby and Froude numbers respectively). Indeed, having non-zero Fr 
and/or Ro opens the door for near-resonances to enter the dynamics and these interactions can, 
both quantitatively and qualitatively, change the behavior of the system [HI . This naturally mo- 
tivates the construction of a model (called the GGG model) that includes all possible wave-mode 
interactions. Quite interestingly, as the two-dimensional (2D) stratified Boussinesq system only 
supports wave modes, this new 3D system can be looked upon as an extension of the full 2D 
stratified problem [20] , [2T] , [22] . 

So far, our discussion has been restricted to considering wave modes in isolation. In reality 
rotating and stratified fluids support an additional vortical mode of motion p3], and in fact in 
geophysically relevant limits, interactions among vortical modes lead to the celebrated quasi- 
geostrophic (QG) equations [l3] , [19] , [23] . Indeed, one of the many contributions by Prof. Majda 
(in collaboration with Prof. P. Embid) in the field of geophysical fluid dynamics concerns a formal 
statement on the emergence of QG dynamics, from the governing rotating Boussinesq equations, 
in the limit Ro ~ Fr = e ^ while holding Bu ~ 1 [25], [19]. In addition to clarifying the 
nature of the QG equations and constructing a general framework for averaging over fast-modes 
in geophysical systems [26] , their work also predicted the emergence of a new regime, the so-called 
vertically sheared horizontal flows (VSHF) when Ro ~ 1 while Fr = e — > j25j.|27j. Since then, 
numerical work in the appropriate parameter regime has confirmed the emergence of VSHF modes 
when considering a rotating Boussinesq fluid under random forcing |24] .128] .1291 .130]. 

As it happens, the vortical mode plays an important role in the redistribution of wave-mode 
energy by means of the so-called wave- vortical- wave interactions. In fact, this is thought to be the 
primary manner in which energy is transferred from large to small scales in rotating and stratified 
flows, and is attributed to be the fundamental mechanism behind geostrophic adjustment |13] . 
Further, in both decaying and forced scenarios, numerical simulations of rapidly rotating and 
strongly stratified (i.e. small Ro, Fr) flows have shown the equilibration of wave modes (by means 
of the aforementioned forward transfer) while energy continues to be transferred upscale in the 
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vortical modes via vortical- vortical- vortical (QG) interactions |13].|14j. For f = N (where /, are 
the Coriolis parameter and the Brunt- Vaisala frequency) in a unit aspect ratio, the wave modes 
are passively driven by the vortical modes. For this case f = N, there are no wave-wave-wave 
resonances (or near-resonances), and the distribution of energy among the wave modes is strongly 
influenced by the presence of the vortical mode [l3],[ll]. Of course, in the atmosphere-ocean 
system f ^ N, leading to the possibility of near and exact wave- wave- wave resonances, though as 
stated earlier these sets of interactions are sparse [19] • All in all, this line of reasoning leads one to 
believe that the vortical mode should in fact play a significant role in determining the wave-mode 
energy distribution. 

These scenarios then lead to a somewhat dichotomous state of affairs with regard to IG wave- 
mode energy, i.e. on one hand we have theories based on the interactions of IG modes in isolation. 
On the other hand, there is reason to attribute a prominent role to the vortical mode in the 
overall dynamics of a rotating and stratified fluid, which includes the distribution of energy among 
the wave modes themselves. This dichotomy — in the context of purely stratified flows — has 
been pointed out and investigated by Waite & Bartello [30]. By forcing only wave modes of 
the Boussinesq system, they inquired into the turbulence generated by a stratified fluid. Their 
principal conclusion was to highlight the importance of the vortical mode in mediating wave-mode 
energy transfer and distribution. Further, they noted the difficulty in properly resolving wave-mode 
interactions [30j. In the present work, in addition to studying the full Boussinesq equations, we 
derive a model that consists of only wave-mode interactions, which in principle provides the means 
for determining the relative importance of wave- vortical- wave and wave- wave- wave interactions for 
the distribution of wave-mode energy. Specifically, by simulating the newly derived GGG model, 
we attempt to see if this system is capable of a robust forward transfer of energy, and compare the 
energy spectra of wave modes from the GGG model to the corresponding wave mode spectra that 
arise from simulations of the full rotating-stratified Boussinesq system. We proceed to introduce 
the basic equations, review the possible means of wave mode evolution and derive the GGG model. 
We then describe the numerical setup and proceed to the simulations of this system. Finally, we 
collect and discuss the results. 

II. BASIC EQUATIONS 

The rotating Boussinesq equations in a 3D periodic setting are [11] 
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^ + fz xu = -Vlp - NOz + z^V^u 

-Nw = kV^9 , V • u = (1) 

where u = {u,v,w) is the 3D velocity field, ip = p/po is a scaled pressure and / is the Coriolis 
parameter with rotation assumed to be aligned with the z-axis. Equations ([T]) result from 
considering periodic perturbations to a state of hydrostatic balance wherein the density profile 
satisfies p = po — bz + p' with <C Po- Further, we have set p' = {Npo/g)0 where 

N = {gb/pQ^/^. Finally, u and K are the viscosity and diffusivity respectively. 



Linearizing (|lj) about a state of rest, substituting solutions of the form = (/)exp{ik-x — cr(k)t} 
one finds 

Here k is the magnitude of the wave- vector and kf^ = k"^ + ky. Modes corresponding to a± are the 
IG waves (referred to as wave modes), while the mode corresponding to gq is referred to as the 
vortical (or geostrophic) mode |23].|13j.|24j. Denoting the eigenfunctions corresponding to fTo,c7± 
by , (j)± respectively, it is known that these eigenfunctions are mutually orthogonal and form a 
complete basis (see Embid & Majda [26], j25j and Smith & WalefFe [23j for details). This enables us 
to project a Fourier expansion of the solution to ([T]) onto this eigenfunction basis, i.e. projections 
of the form 

v(k, t) = ao(k, t)0o(k) + a+(k, t)0+(k) + a_(k, t)<^_(k), (3) 
where v is the state vector, v-^ = {u,v,w,6) |25j . |13] . |23j . Symbolically this yields 

+ ia„(k)a«(k, t) = Y,Yl cS>iP, t)aM, t) (4) 

A /3,7 

where A represents a sum over p, q such that k = p + q. The indices a, (3, 7 run over 0, +, — and 
^kpq interaction coefficient. The exact form of the interaction coefficients is not emphasized 

here, but it is helpful in our discussion to think in terms of the following symmetric definition 
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where contains the first three components of (jijs corresponding to the velocity entries and 
(•)* denotes complex-conjugate. Notice that, with this symmetric definition, C^pq = C^qp! 
coefficient C^pq involves the advection of the mode corresponding to q by the mode corresponding 
to p and vice versa. In the following, a triplet (•, •, •) refers to a particular type of triad interaction, 
and a v refers to a vortical mode while a w refers to a wave mode (either a + or a — mode). For 
example, {w, w, v) (with all permutations implied) represents all possible interactions between two 
wave modes and one vortical mode. 

A. Wave-mode evolution 

An alternate way to represent the expansion is v(k, t) = 6o(k, t)(/'o(k) exp[— icJo(k)i] + 
6+(k, t)(/)+(k) exp[— i(7+(k)t] + 6_(k, t)</)_(k) exp[— i(j_(k)t]. In this case (|4| takes the form 



E E CkS^/3(P, *)^(q, t) exp{i[a(k) - a{p) - aid)] t} (6) 



A /3,7 

which shows the influence of the dispersion on the advective nonlinearity. Indeed, it is now possible 
to construct a dynamical hierarchy consisting of (i) resonances, i.e. |o"(k) — cr{p) — cr{q)\ = 0, (ii) 
near- resonances, i.e. < |cr(k) — 0"(p) — 0"(q)| < 6 where 6 = min {Ro, Fr) and (iii) non-resonant 
interactions, i.e. |iT(k) — cr(p) — cr(c[)\ > 5. In this hierarchy it is postulated that resonances are 
always important, near-resonances play a role when < 6 < 1 i.e. rotation and stratification are 
strong, but we are outside limiting scenarios where 6^0. Finally, for rapid rotation and strong 
stratification, non-resonant interactions are only expected to contribute at higher-orders (in 5). 
While it is not possible for non-resonant interactions to lead to qualitatively different behavior as 
compared to near-resonances, we note that near-resonances open up pathways that are precluded 
under resonance (for details see Sukhatme & Smith [Ij] ) . An explicit example of such a qualitative 
change induced by near-resonant interactions is provided via cyclonic-anticyclonic asymmetry in 
rotating flows [31]. Other examples, such as the influence on the statistical features of water 
waves and capillary wave generation may be found in Annenkov &: Shrira [32j, and Watson & 
Buschbaum [SSj, respectively. 

Considering only resonant and near-resonant interactions, the wave modes can only evolve via 
the {w, w, v) and (w, w, w) classes (see Bartello [E] or Sukhatme & Smith [14j for a detailed view of 
wave and vortical mode evolution). When 1/2 < f /N < 2 the {w, w, w) (near) resonant interactions 
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are absent and — within this dynamical hierarchy — the distribution of wave-mode energy is 
determined by {w,w,v), indicating the crucial role played by the vortical mode. Further, it was 
shown by Bartello [13], that purely resonant {w, w, v) interactions are catalytic in that they involve 
the exchange of energy among two wave modes by using the vortical mode as a catalyst. This leads 
to the elegant result that when f = N the wave modes are passively driven by the vortical mode 
|13| (this follows as near-resonances are not possible for a rapidly rotating and strongly stratified 
fluid when f = N [H]). Of course, in the most general case we lift the restriction 1/2 < f /N < 2. 
As both {w,w,v) and {w,w,w) are viable means of affecting the wave-mode energy distribution, 
it is natural to inquire into their nature and importance. Specifically, given the formal sparsity 
of {w,w,w) interactions [19], are (w^w^v) always dominant? Indeed, do the {w,w,w) interactions 
play any role in the 3D rotating Boussinesq system? Also, what is the nature of these interactions, 
i.e. do they yield a forward transfer of energy under large-scale forcing? If so, do they lead to 
an equilibration of wave modes and yield a well-defined power spectrum? To inquire into these 
questions, we construct a model (the GGG model) that includes only wave-mode interactions. In 
principle, comparing the results from the GGG model with the full rotating Boussinesq system 
should help in settling some of these issues. 

B. The GGG model 

We now present the GGG model in physical space and its derivation. It is useful to first 
determine the mode amplitudes in terms of the physical variables. Towards this goal, begin by 
defining the stream function ip, the potential function x, and the horizontally averaged flows u{z) 
and v(z) in the standard way 



u = Xx - i^y + u{z) and v = Xy + i^y + v{z)- (7) 

Here (•) denotes a horizontal average. These definitions imply V|x = Ux + Vy and V^V = Vx — Uy. 
Substituting ([7]) for u and v in the state vector ,v, ([3| and the orthonormality of the eigenfunctions 
yields (please see the Appendix for additonal details) 
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o+ = — !— (-^5^-/Mrf + iJVM) (9) 

a-, = ^=^(-^^ + fkM - iNk,9) (10) 

< = ^(^)sgn(/fc,) a+ = '^"^ a^^ = (11) 

Next, notice that by adding and and also subtracting from the physical variables 
decouple further, 



< + Ofc = 1 (12) 

^^(y^? + iA:.V^) (13) 
(«£ + «fc J = ^^2^ i< - «L ) = -i V2^- (14) 



In fact, it is now clear from mh and (13) that defining 



M.V?.,-I| and (15, 



is advantageous since this means 



"'^ " < + "fe = ^^fe -^fc- (16) 



In addition M and i? have physical significance : M is the linear potential vorticity (PV) and R 
is a measure of the geostrophic imbalance. It is apparent that the vortical modes contain all the 
linear PV and that the variable R is associated with the IG modes [13,- In the following it will 
prove useful to make the following definitions 



(V2 + ^a,,) ; A = 0-iM, 5 = 0-^R. (17) 



Using these definitions note that (15) can be re-arranged in the form 



(18) 
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A systematic method to complete the derivation of the GGG model in physical space involves the 
three equations obtained from Q, one for each a (0,+,—), after retaining only the interactions 
among the IG modes. Therefore, when a is there are no interactions and the right hand side is 
zero; whereas, when a is + or — , the right hand side is modified by restricting the sum involving 
/3 and 7 to run over + and — only. Then in order to substitute in the physical space variables 
identified with it (above), add and subtract the equations corresponding to a being + 
and — . Finally substitute the physical variables for the mode amplitudes and inverse Fourrier 
transform the equations. This process involves a fair amount of algebra and relevant details are 
presented in the Appendix. 



Here, for the sake of brevity, we present a more illuminating and much shorter derivation. This 
is possible now that the appropriate physical space variables corresponding to the vortical modes 
as well as the IG modes have been identified. It is convenient to rewrite the original system given 
by ([T| in a form that has a separate equation for the time derivative of each of the variables 
corresponding to the modal amplitudes. A consistent system of this sort contains a time derivative 



for M corresponding to the vortical modes and time derivatives for f\/f^R, V w, u{z), and v{z) 



which all correspond to the IG modes (see (12 (13) and (|14[)). This reformulation of (Hj results in 



OMR 
at 



dM f 

— + z-Vx(u.Vu)-^9,[(u.V)0]=O, 

N'^0w + NVl[{u • V)e] + fd,{z • V X (u • Vu)) = 0, 



dV^w 
dt 



+ fVlR + Vliu ■ Vw) - d,{Vh ■ (u • Vu;,)) = 0, 



du{z) 
dt 



fv{z) + dz{uw) = 0, 



dv{z) 
dt 



+ fu{z) + dz{vw) = 0. 



(19) 



Here, the M equation in (19) is obtained by subtracting f /N times the z derivative of the 9 
equation. 



39, 
dt 



+ d,{{u ■ v)e) - Nw, = 0, 



(20) 



from the vertical component of vorticity equation. 



+ Z-V x{u-Vu)-fw, = 0. 



(21) 
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The R equation follows by adding times the the horizontal Laplacian of the 9 equation, 



dt 



+ Vl{{u-V)e)-NVlw = 0, 



(22) 



to / times the z-derivative of the vertical vorticity equation (21). The equation for the Laplacian 



of w in (19 1 involves first solving for if in M by taking the divergence of the momentum equations 



and applying the incompressibility constraint, this yields the formal relation 



^ = V'\-Ne,-V-{u-Vu) + fVl^P). (23) 
Now, taking the Laplacian of the w equation in ([T]) using (23) for (f> gives the relevant equation in 



(19). Note that in taking the required derivatives we lost information about the time derivatives 
of horizontally averaged u and v fields; therefore, we must include the horizontally averaged 
time derivative equations for u and f in ([T| to complete the reformulation (19). Note that the 



horizontally averaged information for 9 is contained in the M equation; the vortical modes for 
kh = do not need separate consideration |37] . 

It is clear that the GGG model, i.e. a dynamical system that only involves wave-mode interac- 
tions should not admit any interactions involving M. This is most easily achieved by using new 
variables {u' ,v' ,w' ,9') where 

2 f2 



^' = Xx- J^S^y + U{z) , v' = Xy + JpSzx + V{z) 



w'^w,9'^ IvlS (24) 



On comparing the primed variables in (24) to the original ones in (19) we see that the new variables 
achieve the goal of removing terms involving the vortical mode from the advective nonlinearity. 
Further, since all interactions involving a vortical mode are to be removed, the evolution of M 
itself must be set to zero. This results in the following concise form for the GGG model 
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dM 

ir = » 

N'^0w + NVl[{u • V)^'] + fd,{z • V X (u' • Vu') = 



^f^h^ AT2rA„.. , ATV72 



dt 

«^ + fVlR + VUu' ■ Vw) - 9,(V, • (u' • Vu',)) = 



du{z) 
dt 



fv{z) + d:,{u'w) = 



dv{z) 
dt 



+ fu{z) + d,{v'w) = (25) 



Note that (25) are well-defined PDEs in physical space, in fact they also conserve energy |S3] (see 



Remmel & Smith [35] for a detailed discussion). 



III. NUMERICAL RESULTS & DISCUSSION 



The Rotating Boussinesq Equations and the GGG model, i.e. (19) and (25) respectively, are 
solved using pseudo-spectral codes in a triply periodic box of dimensions H x L x L, implying H 
{L) is the vertical height (horizontal length) of the domain. The time stepping is done using a third 
order Runge-Kutta scheme. Energy is removed at small scales via an eighth-order hyperviscosity 
term, and at large scales via a linear relaxation term. The combination has the following generic 
form in Fourier Space 

^ + = {-^k^'-a)m, (26) 

with a = 0.1 * I{k < 2.5} and 



E{km,t)\ ^ ,2-16 



. = 2.5 ( k-i-^\ (27) 

Km, / 



In (27) kjn is the highest available wavenumber and E{km,t) is the value of the energy spectrum 
at wavenumber km- All linear terms, including the aforementioned diffusion terms, were treated 
with an integrating factor. 

The numerical results were obtained by forcing each system from a state of rest. The forcing is 
uncorrelated in time, with wavenumber spectrum F(k) given by a Gaussian with standard deviation 
s = 0.5 and energy input rate e = 1, 

v27rs 
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The peak wavenumber kf is chosen as kf = 4 for all results reported herein. The non-dimensional 
parameters are defined as Fr = U/{NLz) and Ro = U/{fLh), where U is derived from the 
energy input rate and the scale at which the energy enters the system, i.e. U = (e/kj)^/'^. 
The length scales {Lz,Lh) are consistent with U, i.e. = ^/kf and Lh follows from the 
prescribed aspect ratio Lh = L * Lz/H. As in atmosphere-ocean phenomena, we set H/L < 1 
and f < N. In fact, we choose f,N to offset the skewed aspect ratio {H/L) so as to satisfy 
Fr = Ro Bu = Ro^ jFr^ = 1. The numerical resolution of the simulations is chosen to yield 
an isotropic grid, for example if HjL = 1/5, we use five times as many points in the horizontal 
directions. 

The modal decomposition in ([s]) is used to feed energy into the system by directly forcing the 
mode amplitudes. The full system is forced by one of the following two methods : 



1. All three mode amplitudes oq, a± are forced with equal weight (see Q). 

2. The amplitude ao is set to zero at every time step such that only wave modes are excited. 

The latter forcing is always used with the GGG model since there are no interactions involving 



vortical modes. Note that for fixed e = 1 in (28), the computed energy input rate for method 2 



will be reduced by a factor 2/3 as compared to method 1, changing the Ro and Fr numbers by a 
factor of (2/3)^/'^ = 0.87. We present the comparison of GGG and full model simulations keeping 
the same level of forcing for the wave-mode amplitudes, rather than the same Ro and Fr numbers. 
Thus, apart from the actual random numbers, the forcing of the IG modes is always identical for 
the two methods and for the two systems (full and GGG). The results will be interpreted keeping 
this choice in mind. 



Given our present computational resources, H/L = 1/5 is the most skewed aspect ratio we stud- 
ied (estimates for mid-latitude atmosphere-ocean dynamics lead to 1/10 < f/N ~ H/L < 1/100, 
but unfortunately these small ratios are well beyond our computational capabilities). Results 
obtained using H/L = 1/5 are presented first. The full system, forced using method 1, is studied 
at resolutions of 80 x 400 x 400 and 100 x 500 x 500 under rapid rotation and strong stratification 
(specifically, Ro = Fr = 0.05). As is evident in Fig. ([T]), the wave-mode energy saturates after 
about 5 (dimensional) time units, and is then almost flat for the next 25 time units. In this 
parameter regime, one expects the energy to be transferred upscale among the vortical modes. 
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however the large-scale damping allows these modes to equilibrate (see [TU], [13], |14j for analytical 
and numerical work on the nature of the various mode interactions when Bu ~ 0(1))- We choose 
to halt the simulations as soon as the vortical modes near equilibration, because from then onwards 
the damping directly affects the dynamics. Similarly, if one does not include large-scale damping, 
then the build-up of energy in the vortical modes soon leads to finite-size effects. In essence, one 
has a small temporal window (depending on the forcing scale and the Fr, Ro numbers) during 
which to study the characteristics of the different modes interactions in the absence of either 
large-scale damping or finite-size effects. Focussing on the wave modes, we see that the energy 
transfer is downscale. Further, on equilibration the wave modes yield a well-defined smooth 
power-law. Given our modest resolution, we see in Fig. ([2]) that the wave-mode spectrum is 
consistent with a power law between to k~^^^ for kf < k < k^, where kf,kd are the forcing 
and dissipation wavenumbers. These results are consistent with previous forced and decaying 
simulations with Bu ~ 1 [l3] , [l3| , [36] . Also, note that the results appear to have converged, i.e. 
we do not see any dependence of the scaling and behavior of the modes for different numerical 
resolutions. Fig. ^ shows the wave, vortical and total energy spectra for the two aforementioned 
resolutions. 

It is well-known that the 2D stratified problem (which only supports wave modes), on 
large-scale forcing, yields a clear forward transfer of energy [20],|2T]. Further, the kinetic equations 
derived from considering resonant interactions of wave modes support power-law solutions [T7j , |18] . 
Keeping in mind that the GGG model can be looked upon as a 3D extension of the 2D stratified 
problem, these two pieces of information lead to the possibility that the GGG model will in fact 
support a forward transfer of energy and yield an equilibrated power-law spectrum. Forcing the 
GGG model with method 2, we show the total energy (equal to the wave-mode energy) vs. time 
in Fig. (Q. Quite clearly, the energy does not equilibrate. On comparing the spectra (not shown), 
none of the runs showed a well-defined power-law. In essence, there is a forward transfer of energy 
in the GGG model, but it is highly inefficient compared to the full system forced by method 1. 
This inefficiency manifests itself in the lack of a systematic forward transfer and the resulting 
lack of an (even approximate) inertial range as is observed in the full system. For the given 
set of governing parameters, we speculate that the GGG model is probably under-resolved at a 
resolution when the full rotating Boussinesq system forced by method 1 yields consistent results 
across a range of resolutions. 
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FIG. 1: The wave-mode, vortical-mode and total energy vs. time for the full system. The simulations are 
at resolutions of 80 x 400 x 400 and 100 x 500 x 500 in an aspect ratio of 1/5 and Ro = Fr = 0.05. The 
runs are halted when the vortical modes near equilibration. 

It is interesting to observe how the presence of the interactions involving vortical modes 
compares to the results obtained from the GGG model when the full system is forced identically 
to the GGG model by method 2. Like the energy in the GGG model, the energy in the full system 
does not equilibriate Fig. ([5|. However, the full system does appear do a better job of moving 
the energy to small scales. After an initial growth of energy in the vortical modes (~ t = 5), 
the energy in the full system grows at a slower rate than in the GGG system. Like the GGG 
system, the full system now produces spectra (not shown) that are not well defined power-laws. 
In fact, the picture that emerges is that there is some transfer from the wave to vortical modes 
that results in the activation of {w,v,w) interactions which help the full system move energy 
toward small scales. But, much like the GGG model the {w,w,w) interactions are incapable 
of a robust forward transfer — possibily due to their being under-resolved in this set of parameters. 



The notion that the GGG model supports a forward transfer of energy, albeit an inefficient 
one, but does not yield a steady state, is somewhat unsatisfactory. There is nothing intrinsically 
wrong with such a state of affairs, but given the possibility of their being under-resolved in the 
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FIG. 2: The wave-mode, vortical-mode and total energy spectrum (at t « 28). The results are from the 
100 X 500 X 500 run. Note that the scaling of the vortical modes is consistent with a enstrophy-cascading 
fc"'^ form for kj < k < fc^. The wave-mode spectrum is best described by a fc~*/^ scaling — we do not 
attribute any special significance to this number. Rather, it is drawn to show that the wave-mode scaling is 
consistent with a power law between —1 and —5/3. 

previous simulation, to gain more confidence in the GGG model we repeated the aforementioned 
experiment at higher Fr, Ro (i.e. under weaker rotation and milder stratification) and in a less 
skewed domain. Specifically, setting H/L = 1/3 and Ro = Fr = 0.1, we performed simulations at 
resolutions of 144 x 432 x 432 and 162 x 486 x 486. Note that the aspect ratio 1/3 case allows 
us to use a higher numerical resolution as compared to the aspect ratio 1/5 case. Fig. ([6| shows 
the energy in time for the GGG model forced by method 2 and the full model forced by method 
1. In contrast to Fig. (j4]), now the GGG energy levels-ofF fairly rapidly. Further, as is shown in 
Fig. ([7|, the spectral distribution of the GGG energy is in the form of a smooth power-law. It 
is interesting to note, from Fig. ([7]), that the GGG spectrum is steeper than the full system's 
wave- mode spectrum. In fact, this steepness is qualitatively consistent with the resonant kinetic 
equation solutions [T7] , [18] . 

As the GGG model equilibrates in these milder conditions, it is of interest to explore the 
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FIG. 3: A comparison of the spectra for the full system at the two different resolutions. The simulations 
appear to have converged and are insensitive to the increasing resolution. 

response of the full rotating Boussinesq system forced the same as the GGG model by method 
2. Note that when the full system is forced by method 1 there was a natural flow of energy into 
the vortical mode that allowed the vortical mode to mediate energy transfer between two wave 
modes by means of wave-vortical-wave interactions. Given that wave to vortical transfers are weak 
(i.e. only due to near-resonances), forcing only the wave-modes allows us to see if the {w,w,w) 
interactions in the full system behave in accord with the GGG model. This is precisely the case, 
and along with the attainment of a steady state (Fig. ([S])), both simulations yield almost identical 
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Time 

FIG. 4: The total energy vs. time for the GGG model. The simulations are at resolutions of 80 x 400 x 400 
and 100 X 500 x 500 in an aspect ratio of 1/5 and Ro = Fr = 0.05. Quite clearly, the system is not near 
equilibration; some of the runs were carried out significantly longer times and yet showed no sign of settling 
down. The spectra (not shown) are also evolving in time and do not show smooth power-law scaling. 

power-law wave-mode spectra for Ro = Fr = 0.1 with H/L = 1/3 (see Fig (|9|). 

In essence, this confirms our speculation that that the wave-mode interactions were under- 
resolved in our numerical simulations of rapid rotation and strong stratification (specifically, 
Ro,Fr < 0.05) in the more skewed domain {H/L < 1/5). Quite naturally, the wave-mode 
energy distribution in the full rotating Boussinesq system is then primarily controlled by the 
{■w,w,v) interactions [30,. As also noted by Waite and Bartello [30], our computations suggest 
the possibility that present-day numerical models may be deficient in their representation of 
{w, w, w) interactions when trying to work in geophysically relevant parameter regimes. Further, 
in milder conditions (i.e. Ro,Fr = 0.1 and H/L = 1/3), the GGG model resulted in a steady 
state with a well-defined power-law energy distribution. This leads us to the possibility that, on 
proper resolution, the {w, w, w) interactions may play a role in the wave-mode energy transfer and 
distribution in realistic parameter regimes. Future work with better computational resources will 
inquire into this possibility as well as universality in the GGG spectral scaling. 
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FIG. 5: The total energy vs. time for the full system forced by method 2. The resolution is 100 x 500 x 500 
in an aspect ratio of 1/5 and Ro ~ Fr ~ 0.05. The system does not equilibriate as it did when forced by 
method 1. However, due to the presence of the interactions among the vortical modes it does a better job at 
moving energy to small scales than the identically forced GGG model does. It is possible that both systems 
are under-resolved in this parameter regime. 



IV. CONCLUSIONS 



In rotating and stratified flows, the forward transfer and distribution of energy among wave 
modes is controlled primarily by wave-vortical-wave {w^ w, v) and wave-wave-wave {w, w, w) 
interactions. However, the relative importance of these two classes of interactions is somewhat 
unclear. On one hand, prior analytical [19] and numerical [l3],[ll] work suggests that the 
wave- vortical- wave class of interactions is of primary importance — when 1/2 < f/N < 2, this 
is immediately evident [30]. On the other hand, especially outside the range 1/2 < f/N < 2 (as 
is true in most geophysical scenarios), the kinetic equation approach considering only resonant 
IG wave-mode interactions yields solutions that are consistent with certain observations |17j.|18|. 
In order to inquire into this somewhat dichotomous situation, we constructed a reduced model 
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FIG. 6: Left : 144 x 432 x 432. Right : 162 x 486 x 486. The total, vortical and wave-mode energy vs. time 
for the full model and the GGG model in an aspect ratio of 1/3 with Ro = Fr — 0.1. The behavior of the 
full model is similar to Fig. (jlj. While, in contrast to Fig. Q, now the GGG energy does level-off. Notice 
the agreement for both the full system and GGG between the two resolutions. 




FIG. 7: Left : 144 x 432 x 432 (w t = 17.5s). Right : 162 x 486 x 486 (« t = 15.5s). Wave-mode energy 
spectra for full and GGG models in an aspect ratio 1 /3 with Ro = Fr = 0.1. These are spectra from a single 
time snapshot of the models in the quasi-statistically steady range. Quite clearly, both may be interpreted 
as fairly clear power-laws with the GGG model spectrum being considerably steeper than the full system's 
wave-mode spectrum which scales between —1 and —5/3. 



20 



3.5| 1 1 1 . 1 3.5 




Time Time 



FIG. 8: Left : 144 x 432 x 432. Right : 162 x 486 x 486. The GGG and fuU system wave, vortical and 
total energy in an aspect ratio 1/3 and Ro = Fr = 0.1 setting in which identical forcing by the IG modes 
was used. Quite clearly, both systems attain an energetically steady state. Notice the agreement for both 
the full system and GGG between the two resolutions. 




log(k} logfk} 

FIG. 9: Left : 144 x 432 x 432 (w t = 18.65s). Right : 162 x 486 x 486 (w t = 17.7s). The GGG and full 
system wave- mode energy spectrum in an aspect ratio 1/3 and Ro = Fr = 0.1 setting in which identical 
forcing by the IG modes was used. These are spectra from a single time snapshot of the models. Quite 
clearly, both may be interpreted as fairly consistent power-laws. 
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consisting only of wave- mode interactions (resonant and non-resonant). This GGG model is given 
by a PDE that conserves energy in the inviscid, unforced case. We studied the GGG model in an 
idealized periodic setting, though its applicability is not limited to such idealized domains. 

We compared the GGG model to the full system in a geophysically relevant parameter regime. 
Specifically, with a skewed aspect ratio H/L < 1, we chose f^N such that Ro = Fr =^ Bu = 1. 
Further, the flow was rapidly rotating and strongly stratified, i.e. Fr = Ro ^ 1. In the full system, 
as anticipated we observed an inverse transfer of energy among the vortical modes accompanying 
by a power-law that was consistent with a scaling for kf < k < k^. The wave- mode energy was 
transferred to small scales, and these modes equilibrated quite rapidly. In accord with prior work 
|13j.|14j.|36]. the wave-mode energy distribution was consistent with a power law with scaling 
between a k~^ to k~^^'^ form for kf < k < kd- 

For the GGG model, we also anticipated a statistically steady wave-mode energy spectrum 
with power-law scaling. This expectation follows from previous results establishing that (i) power- 
law solutions are supported by the kinetic equations governing resonant wave-mode interactions, 
and (ii) the GGG model is a 3D extension of the 2D stratified problem (known to support a 
robust forward transfer of energy). A clear power law was indeed observed for "mild" conditions 
with H/L = 1/3 and Fr = Ro ~ 0.1. Consistent with the kinetic equations for resonant wave 
interactions, the GGG energy spectra were observed to be steeper than the wave- mode energy 
spectra associated with the full equations when all modes were forced with equal weight. When 
the forcing was restricted to excite only wave modes, then the GGG and full systems yielded 
wave-mode spectra with essentially identical power-law scaling. The results indicate that {w, w, w) 
interactions can play a significant role in the transfer of energy from forced scales to smaller scales, 
especially when the wave modes are preferentially excited. However, the energy did not equilibrate 
for the GGG flow in the smaller aspect ratio H/L = 1/5 with stronger rotation and stratification 
{Fr = Ro = 0.05), suggesting an in-efficient forward transfer of energy. Furthermore, the full 
system also did not equilibrate for these parameter values when only the wave modes were excited 
at large scales. These results serve as a caution that improper resolution of wave-mode interactions 
may be a significant issue in present-day numerical models that attempt to work in geophysically 
relevant parameter regimes, consistent with |3D]. Finally, there are a number of issues that raised 
herein that we feel merit further examination. For example, is there universality in the wave-mode 
spectrum from the GGG model (with respect to decreasing Fr = Ro in a fixed aspect ratio). 
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Similarly, with better numerical resomxes, an important issue to be addressed is the difference 
in scaling of the wave-mode energy from full model simulations (with forcing of type 1) and that 
from the GGG model. This gets to the heart of the matter with regard to energy distribution 
resulting from {w, v, w) and {w^ w, w) interactions respectively — indeed, the intriguing possibility 
that the scaling from these two scenarios may converge with progressively stronger rotation and 
stratification provides adequate motivation for such an exploration. 
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V. APPENDIX 



Here we present the GGG model by utilizing the eigenfunction decomposition suggested in Q. 
For the general case, i.e. k^, ky,kz ^ 0, the eigenfunctions are 
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here (•)* denotes complex-conjugate. The special case, kh = 0, kz ^ (the VSHF mode) is treated 
separately. Here = 0, o"='=(k) = it/ and 
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Substituting in ([3| and using the orthogonality of the modes, with as the streamfunction, we 
have 



For kr,,ky,kz 7^ : 



1 



V2akk kh 

_l^_aj^ + fk.kn^l^ - \NkJ) 
V^CTkk kh 

, _ V2kwk 



^^(y^ + ifc.V') (32) 



For kh = 0,k,^0: 



zrr ^ + -1''^ + ^^ - iu + v 

«fc, = 0[z)sgn{fk,) al = ^ a^^ = — ^ 

« + «fc J = («t - «fc J = (33) 



where (•) denotes a horizontal average. 



Using (29) and (32), the GGG model (retaining only IG mode interactions) for the kx,ky,kz / 
case is 



-MO-kl«k "k^ - Z^k=p+q yaafefc I ^ ^pP n -^-jg J\ Phqukh 
PhQhkh PhQhkh PhQhkh ' 
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and 
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Phqhkh / \ "pP ) \ "il / I PhQhkh Ph<lhkh 

ip(^pPzqUPh-<ih)ql , if^^pq'zP'iql , '^Pf^pPzqlivh-^ihip'i , '^popqi{vh-<\h)p\ 

Phqhkh Phqtikh Phqhkh Phqhkh 

The VSHF contributions to the above are 



dt 



+« - VJ(^)(-9-^'' + i.q'. + 9x.gD} (36) 

and 



d{a+ - Ok 



and finahy, using (30) and (34), the evolution the VSHF modes themselves is given by 



+('^)("N=?^)(/P.P.^. + /Pxg.2)| (38) 



and 



H'^) C^''-''T''"''^ ){-^Pyf{p'z+Pzqz)) 
+("^)( "^"".;;^""" )ip./(g.^+P.g.)|- (39) 
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It can be verified that inverse transforming these equations leads to exactly the set specified in 
Section IIB. 
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